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Thermal Modeling of a Single 
Ultralight Fabric Reflux Tube 


Executive Summary 


A thermal model of a single Ultralight Fabric Reflux Tube (UFRT) was 
constructed and tested against data for an array of such tubes tested in the 
NASA-JSC facility. Modifications to the single fin model were necessary to 
accommodate the change in radiation shape factors due to adjacent tubes. There 
was good agreement between the test data and data generated for the same cases 
bv the thermal model. The thermal model was also used to generate single and 
linear array data for the lunar environment (the primary difference between the test 
and lunar data was due to lunar gravity). The model was also used to optimize the 
linear spacing of the reflux tubes in an array. The optimal spacing of the tubes was 
recommended to be about 5 tube diameters based on maximizing the heat transfer 
per unit mass. The model also showed that the thermal conductivity of the Nextel 
fabric was the major limitation to the heat transfer. This led to a suggestion that the 
feasibility of jacketing the Nextel fiber bundles with copper strands be investigated. 
This jacketing arrangement was estimated to be able to double ' tb ® 0 ! he ^ ma ... „ .. . 
conductivity of the fabric at a volume concentration of about 12-14 /<>. Doubling the 
thermal conductivity of the fabric would double the amount of heat transferred at 
the same steam saturation temperature. 


The thermal model consisted of a one-dimensional, composite 
convection-conduction model. The one-dimensional assumption was very good 
on the shaded side of the tube that experienced condensation and was less so 
when the other side experienced solar insolation of sufficient magnitude to cause 
heat addition to the tube. The radiation was modeled as involving diffuse surfaces. 
On the shade side radiation exchange with the lunar surface and space was 
considered. On the sun side, a solar heat flux was input at the radiosity node of 
the fin since it really represented the net energy exchange with the sun. 


Introduction 

A project was initiated at Lamar University to design and test a thermal 
model of the new type heat transfer fin shown in the sketch below. The model 
results were compared with data obtained in NASA’s facility at JSC (Building 33, 
Chamber E) and then used to predict performance on the lunar surface. 
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Nextel 0.56 mm thick 
Titanium 0.0026 in thick 


Condenser 

Section A-A j Section 

I I 1.13 m 

A A 


| 5.1 cm | 

Evaporator 

Section 

Figure 1. A sketch of the Ultralight Fabric Reflux Tube (UFRT) to be modeled. 
The Model 

The thermal model consisted of a sequence of quasi-one-dimensional 
processes: starting at the saturation temperature of the steam inside the tube the 
processes are; a convection process at the tube inside wall, a conduction process 
through the titanium liner, a conduction process through the containment fabric, 
and a radiation process from the fabric outer surface. The process is shown by 
the Sketch of the electric circuit analogy shown in Figure 2. below. 





°-^VvVv -0 “^MA' -0 "^NMAr" 0,— — ~ 

Rcond ^ tl ^ Nex 



Figure 2. A sketch of the electric circuit analogy of the thermal model. 
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Rcond represents the condensation process inside the tube. In Appendix ] it was 
shown that the condensation process could be modeled as laminar film 
condensation on a flat plate. Further it was shown that the magnitude of the 
temperature drop across the film for the range of interest of our study was less 
than 0.1 degree Celsius. The heat transfer coefficient for the condensation 
process was taken at a constant value for the JSC test and for the lunar 
environment at the values shown below. 


h JSC =>60 


kW 
m K 


** lunar® 2 


kVV 
m 2 K 


Conservative Estimates from Appendix I 


Unfortunately, it was postulated in Appendix II and demonstrated in analyzing test 

data in Appendix III that a simple one-dimensional model would not adequately 
represent the behavior of the system. Since solar insolation impacts only one side 
of the tube, the model had to be divided into two parts as shown in the sketch of 

Figure 3. below. 



Figure 3. The two parts (sides) of the thermal model. 
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Some combinations of the radiant boundary conditions cause the direction of heat 
transfer to reverse. Such reversal changes the internal convective process from 
condensation to free convection. In Appendix VIII the free convection heat transfer 
coefficients are calculated to be: 


li 


JSC 


= 57 


watt 
m 2 K 


Conservative Estimates from Appendix VIII 


,, wait 

’lunar® \ 

m K 


Comparison with Experiment 

In Appendix III the thermal model was used to generate data for the cases 
of the NASA-JSC test. 


Conclusions 


1. The thermal model contained in this report closely reproduces the test data 
generated in the NASA-JSC test facility. 


2. From a heat transfer perspective, the optimal tube length to diameter ratio does 
not appear to have been reached. It appears that the optimal value would be 
limited by the thickness of the laminar condensate layer at the bottom of the tube 
which could be increased above the values calculated. 


3. A simple model of the tube, heat transfer fluid system, and support structure for 
the fins indicates that the optimal pitch of a linear tube array should be about 5 
tube diameters. This result is based on a mass of heat transfer fluid, fluid lines and 
support structure which would be about three times that of the tube itself per meter 
of pitch spacing. The parameter optimized, heat transfer per Kg of mass, is very 
sensitive to smaller values of pitch to diameter ratio but much less sensitive to 
larger values (see Appendix VI ). 


4 As expected, the thermal conductivity of the fabric containment limits the heat 
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transfer from the tube. A method of introducing a small volume fraction (about 
12-14%) of copper fibers into the Nextel fabric is suggested in APP^ 1 ^ which 
it is estimatedwould double the effective thermal conductivity of the fabric and 
thereby double the heat transfer from the tube at the same saturation tempera 
of the steam inside the tubes. 



End of Section 
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APPENDIX I 


Thermal Modeling of the Condensation Process 
Inside a Single Ultralight Fabric Reflux Tube (UFRT) 
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Project Objective: To formulate and test a thermal model of the new type heat 
transfer fin shown in the sketch below. The model is to be tested against data 
obtained in NASA's facility at JSC and then used to predict performance on the 

Lunar surface. 



Section A-A 


Nextel 0 56 mm thick 
Titanium 0.0026 in thick 
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Ti-2AI-2Mn 


Vol 4 20 cm 

VWI water 


1.13 m 


The first step will be to examine the condensation process in the condenser 
to determine the complexity of the model needed. We will begin by making 
simple assumptions of laminar filmwise condensation on a vertical flat plate 
with no interface shear between the vapor and liquid film. We must test to 
see whether each of these assumptions is valid. 

To begin we will take saturated steam at 20 Ib/im abs condensing on the 
inside of ’a vertical 1 in diameter tube 1.13 m tall. The tube wall temperature 
is constant and varied from 227.8 to 227.95 <>F. The heat transfer 
coefficient, the heat transfer and the condensate flow will be calculated for 

each wall temperature. Take T sat = 227.96 °F and h fg = 960.1 Btu/lb. 

Ill 

j 0 15 h 227.8 . j .01 g 9.8 ^ 

sec 
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Heat Transfer - Watts 





Bottom Laver Thickness - 







Bottom Laver Thickness - 


0.045 





Revnold's Number 
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Therefore, we can examine our assumptions. The filmwise condensation 
process is laminar as shown by Reynold's Numbers in the range of 2 to 14 
in the figure above. The figure for condensate layer thickness at the 
bottom of the tube gives a maximum of: 

^ max 3.5 10 ' in 

at the maximum heat dissipation that we are interested in of 150 Watts. 
The ratio of diameter to condensate thickness is then: 

ratio 5 ' naX ratio = 1.378 10 3 

(1 

which validates our vertical flat plate assumption. 

We will calculate and plot the upward vapor velocity to determine 
whether our no interface shear assumption is approximately correct. 

We know that for steady state the upward mass flowrate of vapor must 
be equal to the downward flow of condensate. 

* •' 26 i ) 2 . v "j 

A v. ' . w vapor w condensate v v. . 

.1 4 1 I P Vj A Vj 


Upward Vapor Velocity 
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which validates our assumption of no interface shear. 

We know that the initial water charge is 20 cm 3 . Let us check to see the 
maximum volume of water contained in the film to be sure that the 
evaporator section does not go dry. The volume of water contained in the 
film is given by: 
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2 

7 E d L i 5.1 cm ^ 

vapor 4 m vapor = 2,421 ‘ ,0 kg 

Volume of Water in Film 



So about 10% of the 20 cm 3 charge is contained in the film at the maximum 
anticipated heat load of 150 Watts, however, the mass of water vapor becomes 
the same order of magnitude for a saturation temperature of 500°K. 


Lunar Environment 


Now we need to repeat the calculations for the lunar environment. 

9.8 ni 

g 6 ~2 
” sec 

Again, we will take saturated steam at 20 lb/im abs condensing on the 
inside of a vertical 1 in diameter tube 1.13 m tall. The tube wall temperature 
is constant and varied from 227.8 to 227.95 <>F. The heat transfer 
coefficient, the heat transfer and the condensate flow will be calculated for 

each wall temperature. Take T sat = 227.96 °F and h fg = 960.1 Btu/lb. 
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Therefore, we can examine our assumptions. The filmwise condensation 
process is laminar as shown by Reynold's Numbers in the range of 2 to 10 
in the figure above. The figure for condensate layer thickness at the 
bottom of the tube gives a maximum of: 


at the maximum heat dissipation that we are interested in of 150 Watts. 
The ratio of diameter to condensate thickness is then: 


ratio maX ratio = 1.654*10 J 

d 

which validates our vertical flat plate assumption. 


We will calculate and plot the upward vapor velocity to determine 
whether our no interface shear assumption is approximately correct. 
We know that for steady state the upward mass flowrate of vapor must 
be equal to the downward flow of condensate. 
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Delta Temperature - degrees K 


which validates our assumption of no interface shear. 

We know that the initial water charge is 20 cm 3 . Let us check to see the 
maximum volume of water contained in the film to be sure that the 
evaporator section does not go dry. The volume of water contained in the 
film is given by: 
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So about 25% of the 20 cm 3 charge is contained in the film at the maximum 
anticipated heat load of 150 Watts, however, the mass of water vapor increases 

the total to about 35% at a saturation temperature of 500°K. 
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APPENDIX II 


Radiation Modeling 


II - 1 


10/17/96 



Radiation Heat Transfer 


If solar insolation is present, we might note that the insolation is incident on 
only half the fin area as shown in the sketch below. 

Solar Insolation 


Sun Side 


Shade Side 

The area of the fin will be composed of equal areas on the sun side and 
on the shade side. 

Radiant energy exchange will be assumed to be diffuse. The electric 
circuit analogy for the radiant energy exchange between the fin, the 
lunar surface, the sun and black space can be represented by the 
partial sketch below. 

^bsun 

'I 

I 1 

I 

F F 

■— bfi n ‘-bmoon 

* - yw / <j>- Ms — ■■ ■ / , M;~ 

^ ’moon 

' L 1 A 

I] moon moon 

t 

1 - E 

■ pace 

H space A wee 

II E 

^-bspace 

The reason we have presented only a partial sketch is that we will show 
that it is all that is necessary to determine the net radiant energy 
exchange with the fin. In each of the surface resistances there is an area 
term in the denominator. All other tt^mi will be expressed with the 6 





term in the denominator. All other terms will be expressed with the area of 
the fin in the denominator. Since A space , A moon , and A sun are much much 
larger than A fin , their associated resistances can be taken as zero. In 
addition, we know the net heat transfer with the sun to be the solar 
insolation value. Hence, the modified electric circuit analogy would look 
like that below. 


i - 


,: nn 
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fin 
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fin 
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fin 


Solar Insolation 


^bfin 
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ii 
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i 

Afin ^ Tin moon 

1 

Afin ^ fin-space 


For which we can sum the heat flows into the node J fln and set them to zero. 

bfin fi 
I 1 e fiii 
\ 6 Tin ^ fin 


which applies on the sun side and a similar equation would apply on the 
shade side. 
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I £ fin 
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fin2 ^ bspacc fin2 
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^ bmoon ^ fin2 


\ 


= 0 
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In these equations we know that E bspace =0, E bfjn =aT fjn 4 , E bmoon =aT moon 4 and 
e fln =0.89. The configuration factors may be determined in the following 
manner: consider the sketch below, 


space 



lunar surface 


From the figure above it is obvious that F dA _ moon and F dA . space are both 0.5 
for any position on the surface of the fin (assumes lunar surface is an 
infinite flat plate). Therefore F fin moon and F fin . sp ace are 0.5 for either the 
whole fin area or half the fin area (divided into sun side and shade side). 
The equations above then reduce to on the sun side taking q sun to be the 
energy incident on the fin per unit area times the perpendicular area: 


J fin = e fin ° T 
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(i T 
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moon 
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*1 sun 
A fin 
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and on the shade side: 


fin2 _E fin ^ fin2 1 ( * 

where A fin is half the total fin area, T fjn is the outside temperature of the 
Nextel exposed to the sun, and T fin2 is the Nextel temperature on the 
shade side. 


E fin ) 


a T 



We can now examine the simple but important case of what minimum heat 
transfer is necessary to keep the fin from freezing-up during the lunar night. 
For this time we will take T moon = 88 °K, T sat = 273°K and no solar insolation. 

d 1 in L 1.13 in 

E fi„ «•» ^ sat 273 K T mooB 88 K 

d d i 2 .00056 m i 2 .0026 in d o = 0.027-m 

We will iterate for the outside Nextel surface temperature such that the heat 
loss by radiation is equal to the heat convected and conducted through the 
cylindrical fin wall. We will take only half the fin area so that we will be 
consistent later on when we must consider sun side and shade side. 
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The iteration process is shown in the sketch below. 
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There is a more compact way to do the calculations below, however, 
this method has been used to try to illustrate the calculation 
methodology. 
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Therefore it takes q (ot 2 q j q (ot = 25.871 watt 
to keep the fin from freezing during the lunar night. 

Increasing the saturation temperature to 300°K. T sat 300 K 
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at T sat = 350 K 
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End of Section 
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APPENDIX III 

Thermal Model of the JSC Test 
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JSC test Model 
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Therefore the energy balance at the radiosity node is: 
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Then an energy balance setting the energy transfer through the composite wall 
equal to the energy radiated from the outer surface gives: 
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Modeling of a Single Fin on the Lunar Surface 
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Single Fin Model for the Lunar Environment 
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Test Series 1-5 


To model the different test series, set the model input 
parameters as shown - results plotted below. 
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APPENDIX V 


The Shape Factor Between Adjacent 
Tubes in a Linear Array 
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To determine the shape factor between adjacent tubes in a linear array 
we will treat them as infinitely long. Then the shape factor is simply the 
difference in the sum of the crossed strings minus the uncrossed strings 
divided by the base length. The sketch below shows two tubes of 
diameter, d, in such an array separated by a pitch distance p. 
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Pitch Optimization of Linear Tube Arrays 
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For a Saturation Temperature of 373°K, Lunar 
Temperature of 80°K, and a Solar Angle of 45° 


lloal Transfcr/Mflss vs Pitch/Diamcter 



Pilch/Diamoter - Unoar Spacing 

mass frac=l 
mass IVac=2 
mass ffrac=3 

This figure indicates that a resonable design pitch/diameter ratio would be 
about 5. A higher saturation temperature would perhaps push this ratio down 
to 4 to 4.5 as would more massive structure and piping arrangements. Even 
in these cases, the use of a ratio of 5 would not be costly since the function is 
not as sensitive to the right of the peak as it is to the left. 
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APPENDIX VII 

Equivalent Thermal Conductivity of Concentric Layers 
of Copper and Nextel 
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The model below shows copper wires covering Nextel threads in a composite 
fabric weave. The purpose of the fine copper wires is to increase the thermal 
conductivity of the fabric while minimizing the loss of strength provided by the 
Nextel fibers. The copper layer is shown around the Nextel fibers so that the 
weave of the fabric allows the heat to be conducted along the copper wires or 
strands and thereby carries the heat from the inside surface to the cross layers 
of fibers and from these cross layers to the outside surface as shown in the 
figure below. Copper is used instead of gold to avoid an undesireable decrease 
in outside surface emmisivity. 
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APPENDIX VIII 


Estimation of the Free Convection Heat Transfer 
Coefficient Inside the Vertical Tube 
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Estimation of the free convection heat transfer coefficient between the 
inside tube wall and the saturated steam from the evaporator. This 
condition exists when the radiant boundary conditions produce an 
external surface temperature above the saturated steam temperature. 

The heat transfer coefficient will be taken as that for a vertical flat plate 
even though the tube will not behave in this fashion. This will give a 
larger value which in our case will produce a higher heat transfer to the 
tube and thereby give conservative results. 
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Use h = 57 watt/m 2o K for the JSC test and: 


g 


g 

6 


Gr 


g P AT l 3 



Gr = 7.412- 10 10 


still turbulent 


Nil 


i 

0.15 (Gr Tr) 3 

16 


I 


/ 0.492 \ 
l Pr / 


h 



h - 31.2 


watt 
m 2 K 


Nu = 550.02 


h - 5.495 • 


BTII 
hr ft 2 R 


h = 31 watt/m 2o K for the lunar surface. 
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